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ABSTRACT 

We measure the fraction of Luminous Red Galaxies (LRGs) in dynamically close pairs (-with projected 
separation less than 20 kpc and velocity difference less than 500 km s"^) to estimate the dry merger 
rate for galaxies with -23 < M{r)k+e,z=o.2 + 5 log ft. < -2L5 and 0.45 < z < 0.65 in the 2dF-SDSS LRG 
and QSO (2SLAQ) redshift survey. For galaxies with a luminosity ratio of 1 : 4 or greater we determine a 5(t 
upper limit to the merger fraction of L0% and a merger rate of < 0.8 x 10~^ Mpc^^ Gyr~^ (assuming that 
all pairs merge on the shortest possible timescale set by dynamical friction). This is significantly smaller than 
predicted by theoretical models and suggests that major dry mergers do not contribute to the formation of the 
red sequence at z < 0.7. 

Subject headings: galaxies: interactions — galaxies: formation and evolution 
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1. INTRODUCTION 

In the standard A Cold Dark Matter (ACDM) picture, 
galaxies are assembled via the progressive (hierarchical) 
merger of increasin gly more ma ssive subunits (see, e.g., 
ICole et al.ir2008. : .Neistein & Dekelll200& for recent reviews). 
About 1/2 of the stellar mass in present-day massive {L > L*) 
galaxies is expected to be accreted via major mergers at z < 1 
(e.g.,!be Lucia et al. 2006; De Lucia & Blaizot 2007). Merg- 
ers should therefore be a common occurrence in the life 
of galaxies and have a profound influence on their proper- 
ties (such as star formation, morphology and nuclear activity 
among others). 

Observational and theoretical evidence suggests that an in- 
creasing fraction of mergers at lower redshifts should take 
place between spheroidal or gas-poor galaxies ('dry' merg- 
ers). In their st udy of the morphology of mergin g pairs in a 
CDM universe, i Khochfar & BurkertI (I2003L |2005|) found that 
elliptical galaxies brighter than L* were mainly formed via 
major dry mergers and that gas-rich mergers were only im- 
portant at lower luminosities. Above a 'threshold' mass of 
^ 6.3 X 10^" M0, galaxies are not expected to grow their 
stellar mass by (induced) star formation, but their primary 
mode of mass accret ion at z < 1 is via gas-poor (stellar- 
dominated) mergers (iDe Lucia et al.1 120061: iKlioc hfar & Silk 
I2OO81) . Unhke gas-rich mergers, dry mergers are not expected 
to distur b the observed tight scaling relatio ns for early-type 
galaxies (iBovlan-Kolchin et alJl2005L |2006') and to b etter re- 
produce the internal structure of local ellipticals (Naa b et aP 
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l2006Ll2007h . 

The first examples of dry mergers, in significant numbers, 
were observed as close pairs of galaxies o n the red sequence 
in deep images of high redshift clusters (Ivan Dokkum et all 
1999, 2001), where approximately 50% of galaxies were ex- 
pected to have underg one a dry merger to the present epoch. 
In the local universe, Ivan DokkumI (l2005h estimated that ^ 
30% of z < 0.1 bright elUpticals in flie MUSYC and ND- 
FWS surveys showed residual structural features indicative of 
a dry merger in the 'recent' past. 

iBell et all (l200 6a) used close pairs in the Combo- 17 survey 
having similar photometric redshifts and showing evidence of 
interactions to infer an integrated merger rate of ~ 80% since 
z < 0.8 f or red galaxies with Mb < —20.5 + 5 log h, while 
iLin et al.l (|2b08) used dynamically close pairs in the DEEP2 
survey to derive an integrated dry merger rate of 24% at z < 
1.2 for galaxies w ith — 21 < Mb — 5 log h < —19. However, 
iHsieh et al] (1200 8l) identify close pairs in the RCS survey and 
derive an integrated merger rate of only 6% per Gyr since 
z = 0.8 for gala xies with -25 < - 5 log ft < -20, and 
IWen et alj (|2008) use the same approach as lBell et all (l2006al) 
on LRGs in the SDSS survey to determine a merger rate of 
0.8% p er Gyr at z < 012 for galaxies with Mr < -21.2 + 
5 log h. iBundv et al.l (12009 ) find evidence of few to zero pairs 
of red galaxies at z < 1.2 in the GOODS data. 

The small scale correlation function of LRGs in the SDSS 
at z < 0.36 analyzed by Masje di e t al.. (2006^ 2008) is con- 
sistent with an upper limit of < 1.7% per Gyr to the dry 
merger rate for g alaxies wi t h Mj < —22.75 + 51og/i at 
0.16 < z < 0.30. iBell et all (l2006bh used this method to de- 
rive a merger rate of 4% pe r Gyr for Ma < — 20.5 + 5 log h 
galaxies at 0.4 < z < 0.8. IWhite et al.T(l2007h obtain an in- 
tegrated merger rate of ~ 30% for LRGs a t 0.5 < z < 0.9 
in the NDFWS survey. Finally, IWake et al.1 (120081) apply the 
correlation function method to galaxies the 2SLAQ survey to 
derive a merger rate of 2.4% per Gyr at 0.19 < z < 0.55. 

The galaxy merger rate has been usually measured via the 
fraction of galaxies in close pairs, often with the added re- 
quireme nt of closeness in veloc ity space to cull interlop- 
ers (e.g.. iPatton et al.ll2000L llOOl - hereafter. P00,P02), and 
the fraction of galaxies showing significant asymmet ries in 
their light distribution (e.g., Conselice 2003; Conselice et alj 
l2003h . The rationale behind the former approach is that if a 
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merger is to occur, a companion must be present, and there- 
fore the close pair fraction is related to the future merger rate 
for the galaxies being considered. The latter method relies 
on the observation that if a galaxy has undergone a recent 
merger, it is likely to be morphologically disturbed, and there- 
fore an asymmetric light distribution would be a signpost of a 
recent mer ger. A critical discussion of t hese approaches ma y 
be found in lDe Propris et al.l (l2007h and lGenel et dl (l2008al) . 

In this paper we derive the dynamically close pair frac- 
tion for galaxies in the 2SLAQ survey and determine an up- 
per limit to their merger rate. The benefit of using close 
pairs is that we are able to select major merger candidates 
between galaxies in a specified mass range (by appropriately 
choosing the luminosity of and magnitude difference between 
participating galaxies), derive a merger rate within a speci- 
fied timescale (dependent on the chosen projected and veloc- 
ity separations for the pair members and theoretical simula- 
tions) and identify ongoing merger candidates for later study. 
Galaxy asymmetries tend to be more difficult to interpret in 
this fashion and usually require better quality imaging than 
we have av ailable in our sur vey, especially for dry mergers 
(iBell et alJ l2006a: W en et al.l l2008). to which the 2SLAQ sur- 
vey is most sensitive. Unlike asymmetries, galaxy pairs are 
sensitive to the merger fraction of 'progenitor' halos and this 
quantity may be more directly compared to theoretical models 
dOenel et al. 2008a). 

The structure of this paper is as follows: in the next sec- 
tion we present the data, describe the methodology and derive 
the pair fraction and an upper limit to the dry merger rate at 
the intermediate redshifts sampled by the 2SLAQ survey. We 
then discuss and compare our results in the context of galaxy 
formation models and recent work on the dry merger rate. We 
adopt the latest cosmological parameters with ^Im — 0.27, 
Oa = 0-73 and Hq = 100 km/s/Mpc. Unless otherwise 
stated, all absolute magnitudes quoted in the following are in- 
tended as including a term of +5 log h and all distance mea- 
sures need to be referred to h (to the appropriate power). 

2. THE 2SLAQ SURVEY DATA 

The data used in the 2SLAQ survey consist of LRG s with 
i < 1 9.8 mag. from the original sample of Eisenstein^et_alJ 
(1200 lb . selected by th eir q — r and r — i colors to lie at 
0.45 < 2 < 0.65 (see iFukugita et allll996l for a description 
of the SDSS filter system). Photometry and astrometry for 
the target LRGs are derived from the SDSS Data Release 1 
([York et al. 2000; Abazaiian et al. 2003) with improvements 
from the latest release available at the time of the spectro - 
scopic observations (DR4 - Adelman-McCarthv et al.ll2006h . 
The objects lie in two long strips on the celestial equator, di- 
vided into several disconnected patches, each of which is be- 
tween 10 and 30 deg^ in area, for a total survey coverage of 
182 deg2. 

Spectroscopy for candidate LRGs was carried out at 
the Anglo-Austra lian telescope using the 2dF facility 
dLewis et al.l 120021) . For objects in the main sample (Sam- 
ple 8), it was found that most galaxies are within the speci- 
fied redshift limits and were observed with high spectroscopic 
completeness (typically 87%). A complete description of the 
data can be found in the general survey paper by Cannon et 
al. (2006; hereafter C06). Because we have highly complete 
spectroscopy, we can confirm that at least 90% of our galax- 
ies have K-type or LRG spectra, with no sign of star forma- 
tion, while less than 1% of the sample shows emission lines 
(iRoseboom et al.ll2006h . We can therefore use this sample to 



measure the dry merger rate. 

Following IWake et aTl (l2006h . we computed the absolute 
magnitude for objects with reliable redshifts, k ~\- e corrected 
to the SDSS r band at z = 0.2. This facilitates comparison 
with previous work on galaxy evolution and the merger rate 
from the SDSS (e.g., Masjedi et al. 2006, 2008). Note that no 
correction for internal extinction was applied to these galax- 
ies. 

3. METHODOLOGY 

We calculate close pair statistics following the formalism 
developed by Patton et al. (2000, 2002; hereafter POO, P02) 
for the SSRS2 and CNOC2 surveys. Here, we give an 'algo- 
rithmic' description of the procedures used, and show how we 
apply weights to correct for sources of incompleteness and the 
flux limited nature of the 2SLAQ survey. A fuller description 
of the method can be found in POO, P02. 

3.1. Luminosity of the pair sample 

Let us consider a sample of Ni primary galaxies brighter 
than a limiting absolute magnitude Mi in some volume of 
space, and, in the same volume, a sample of N2 secondary 
galaxies brighter than a limiting absolute magnitude A/2. The 
two samples may coincide (i.e., AIi — M2), as is often the 
case for redshift surveys: we then study 'major' mergers be- 
tween galaxies of approximately similar luminosity. We are 
interested in knowing the fraction of galaxies in the secondary 
sample that are dynamically close to galaxies in the primary 
sample. We define two galaxies to be dynamically close if 
they have a projected separation < 20 /i^^ kpc and a ve- 
locity difference of < 500 km s~^, as used by POO, P02 and 
in subsequent work. 

Of course, pair statistics are usually computed from red- 
shift surveys, which are flux-limited rather than volume lim- 
ited. We therefore need to account for the dependence of pair 
counts on the clustering properties and the mean density of 
galaxies in the sample, and to correct for sources of spatial 
and spectroscopic incom pleteness. Because c lustering is lu- 
minosity dependent (e.g.. lNorberg et ai]|2002l) we follow POO 
and restrict the analysis to a fixed range in luminosity, within 
which clustering properties are not expected to vary signif- 
icantly, by imposing additional bright (Mbright) and faint 
{M faint) limits on the sample. This means that we derive a 
pair fraction (and merger rate) for galaxies within a specified 
range and ratio in luminosity. 

In our case, we select galaxies with 0.45 < z < 0.65 
(where the survey is most complete) and with A'hright — ~23 
and M faint = —21.5 mag. Because the range of luminosities 
we survey is small, we let M faint coincide with M2. We 
also let the primary and secondary samples coincide, to se- 
lect major mergers between galaxies of approximately equal 
luminosities (and make full use of the spectroscopic sample). 
Figure 1 shows the distribution of 2SLAQ galaxies in abso- 
lute magnitude vs. redshift, together with selection Unes in 
magnitude and redshift (see below for details). 

The faint limit allows us to include most 2SLAQ galaxies 
at z > 0.45 and reaches slightly below the luminosity of nor- 
mal L* galaxies at the mean survey redshift (or, similarly, be- 
low the turnover in the 2SLAQ LRG luminosity function of 
IWake et aljr2006l) . We choose the bright limit of Mr = -23 
to avoid the rarer and most luminous galaxies in 2SLAQ, 
which are likely to be more biased, and to mitigate the ef- 
fects of luminosity-dependent clustering (as described above). 
With these choices, we study major mergers, with luminosity 
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Fig. 1. — The distribution of absolute magnitudes for 2SLAQ galaxies 
vs. redshift. The thick solid line is the redshift-dependent absolute magni- 
tude limit for selection Mii^(z) (see equation (1) and related explanation 
in text). The thick dashed lines mark the ideal volume limited box to which 
the observed sample is coirected using the weights described in equations 
(2)-(9). 

ratio > 1 : 4, where the secondary is at least 4 times less 
luminous than the primary galaxy'^, for Mr < —21.5. 

We can translate our luminosity range into a stellar mass 
range using the conversion between r band absolute magni- 
tude and stellar mass by Baldrv et alj ((2006), for typical LRG 
colors (u — r > 3) and wi t h assu mptions as to stellar pop- 
ulations as in iBaldrv et al.l (l2006h . With these choices our 
—23 < Mr < —21.5 luminosity range corresponds to a stel- 
lar mass range of 16-5 x 10^" M©. 

3.2. Density weighting S{z) 

The pair fraction that we actually measure needs to be 
corrected to the value that would be observed in an ideal 
volume-limited survey. In a flux-limited sample, primary 
galaxies at lower redshifts will have a greater likelihood of 
having a secondary companion within the survey than galax- 
ies at higher redshift. We correct for this bias by assigning 
a greater weight to the rarer companions found at the high 
redshift end of the survey. This is carried out by comput- 
ing for each galaxy a weight that renormalizes the sample to 
the density corresponding to a volume limited sample within 
Mbright < M < M2- The weight is calculated by integrating 
the luminosity function over the appropriate ranges in abso- 
lute magnitude. 

At each redshift, we then search for pairs of galaxies within 
Mbright and a redshift-dependent absolute magnitude limit 
Miim{z), which is defined as: 
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-5 log dL{z)-2h-k{z)-e{z)] 
(1) 

where i — 19.8 mag. is the apparent magnitude limit of the 
survey, dL{z) is the luminosity distance and k{z) and e{z) are 
the k and evolutionary corrections. Recall that here M faint is 

Note that of course no sample can be made complete for a given 
luminosity rat io without discarding a large fraction of the data (e.g., 
IPatton & Atfield»20Q8i) , so theoretical comparisons need to take into account 
the observational limitations of this technique. 



set to coincide with M2- The fc + e corrections are taken to 
be the maximal corrections for a galaxy formed at high red- 
shift and undergoing pure passive evolution, as other choices 
would allow galaxies to fall in and out of the sample accord- 
ing to their star formation histories (POO, P02). Note that here 
max means 'the brightest of rather than the (arithmetically) 
larger quantity. Fig. 1 shows this selection limit as applied to 
2SLAQ data. 

Each secondary galaxy is weighted by the inverse of a selec- 
tion function S{z) defined as the ratio of densities in volume- 
limited vs. flux-limited samples: 
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where L(M) = ioO-4( M-Me)^^ ^(^^^ jj^g lrg 
luminosity function from IWake et aH (|2006). The integrals 
run from M^right to either the faint absolute limit M2 (in the 
denominator) or to the redshift dependent absolute magnitude 
limit Miim{z) (in the numerator) set by the apparent magni- 
tude limit of the survey. 

Galaxies in the primary sample at low redshift will also 
have the largest number of observed companions, while pri- 
maries at higher redshifts will have fewer observed compan- 
ions. This effect is corrected in a similar fashion as for galax- 
ies in the secondary sample, by applying the inverse of the 
secondaries' weight to primaries, i.e Sn{z) and Sl{z). 

3.3. Spatial incompleteness Wb, Wy 

We need to account for pairs missed because one of the 
galaxies falls outside of the survey footprint. A potential com- 
panion may lie beyond the survey limits on the sky or be hid- 
den in the 'shadow' of a bright star, where galaxies cannot be 
detected or the spectra are contaminated". For each primary 
galaxy we compute the fraction 1 — /b of the Trr^ area that may 
lie outside of the effective survey area. The weight to be ap- 
plied to the secondary galaxies is then = l//h- Primary 
galaxies may similarly be lost in the survey boundaries and 
the appropriate weight to be applied to the primary sample is 

Wb, ^ fb^ -iWfT/- 

For objects with small separations, the SDSS photometric 
pipeline tends to merge pairs into a single galaxy. For the 
LRG sample studied by Masiedi et al. (2006) the pipeline be- 
comes unreliable for < 3". The 20 h^^ kpc search radius 
we use corresponds to angular separations of 4.91" to 4.07" 
at 0.45 < z < 0.65. We inspected all our 7889 images to ver- 
ify whether the SDSS pipeline correctly identifies photomet- 
ric companions, using the SDSS 'Image List/Navigate' tools. 
The results are that the pipeline certainly finds all targets with 
projected separations between 3" and 5". For smaller sepa- 
rations the pipeline is more erratic, especially for fainter sys- 
tems. We therefore exclude an area of 3" around each target 
from our search for companions and correct for the missed 
region via the Wb weight, where we assume that the distribu- 
tion of pairs is uniform with projected separation over these 

" This is calculated using a formula developed by I. Strateva (unpub- 
lished) for the SDSS survey 
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small scales. The size of the added factor lies between 37% 
and 54% depending on redshift. 

The SDSS pipeline also tends to 'double count' light for 
close pairs, making these galaxies brighter than they would 
otherwise be. Because fainter galaxies have more merg- 
ers an d minor mergers are more common ( Pat ton & Atfieldl 
I2008h. this has the effect of artificially raising the pair frac- 
tion. i Masiedi et al.l (l2006h estimate that this should decrease 
the actual merger rate for galaxies with < 5" by a factor 
of about 5. Since this is an uncertain correction and we are 
interested in an upper limit to the merger rate, we do not con- 
sider this effect here (the effect cannot be well modelled for 
our sample, as the SDSS pipeline cannot be run locally). 

We also correct for possible companions missed by the red- 
shift 'cuts' we apply to the survey. If the primary lies within 
500 km s^^ of the redshift boundaries, we ignore all compan- 
ions between the primary and the borders and apply a weight 
Wy^ = 2 to all other companions in the opposite 'direction'. 
We also need to apply a similar weight to account for poten- 
tial primaries lost in the redshift boundary. This weight is 
the reciprocal of the weight applied to the secondaries, i.e., 
Wy, = 1/2 

3.4. Spectroscopic incompleteness wg 

As all redshift surveys, 2SLAQ is not complete to its flux 
limit. Our success rate is 87% for all galaxies within Sample 
8 of C06 but the spectroscopic incompleteness may be depen- 
dent on the separation between galaxies. Fibers in the 2dF 
positioner cannot be placed closer than about 25" from each 
other in every single configuration. However, this effect is 
compensated by the overlap between individual 2SLAQ tiles 
and by the fact that fiber configurations were generally re- 
touched halfway through each (typically 4 hours) exposure to 
place fibers assigned to galaxies for which a reliable redshift 
had already been obtained on to a nearby target (see C06 for 
a description of the observations). 

In order to correct for this source of bias we need to es- 
timate the relative incompleteness for close pairs over the 
range of separations of interest (corresponding to 20 kpc) 
and compare it with the incompleteness at large separations, 
where fiber interactions are not important. We follow P02 and 
estimate this weight by computing the ratio between the num- 
ber of pairs between galaxies with redshift information (N^z) 
within the redshift range of interest and the number of pairs 
in the input photometric sample (Npp), which is by definition 
complete, as a function of angular separation 6 The weight to 
be applied is the ratio between the spectroscopic and photo- 
metric pair completeness at each separation normalized to the 
value at large separations. 

In 2SLAQ we have higher overall completeness than 
CNOC2 (by about a factor of 2) but sample a smaller range of 
projected separations (because of our higher mean redshift). 
We plot Nzzi0)/Npp{0) vs. d in Fig. 2 for Sample 8 targets; 
the error bars are assumed to be Poissonian. Although the 
data are noisy (because there are relatively few potential pairs 
to start with), the value of N^z/Npp al < 6" is consistent 
with the value of this ratio at > 100", arguing that we are 
not systematically more incomplete at small separations than 
at large ones where we are not affected by fiber collisions. 

P02 model the incompleteness in the CNOC2 survey by fit- 
ting a polynomial to the ratio N^ziS) /Npp (6) as a function of 
projected separation 9. Given the small number statistics and 
noisier nature of our data, it is not fully justified to model the 
completeness as a function of 9 with a polynomial as done in 
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Fig. 2. — The ratio of galaxies in pairs between galaxies with redshifts and 
galaxies in pairs in the photometric sample, as a function of pair separation 
9. 

P02. We therefore adopt an uniform weight of 1 for 2SLAQ 
galaxies over the separation of interest, as we do not appear 
to be systematically more incomplete than at larger separa- 
tions. Although pair completeness drops at 30"-60" separa- 
tions, because of fiber collisions, these large separations are 
not relevant to our study (as pairs with rp > 50 h^^ kpc or 
Av < 1000 km are largely spurious - POO, De Propris et 
al. 2007). 

Given the small number statistics, and the fact we are ulti- 
mately deriving an up per limit, we eventually decided to fol- 
low the approach by iDe Propris et al.l (l2005l) to estimate the 
contribution from close pairs missed because of fiber colli- 
sions. We search around each of our main targets (Sample 8 
galaxies with -23.0 < M (r) < -21.5 and 0.45 < z < 0.65) 
for a companion (lying within rp) in the sample of galaxies 
(from the input photometric sample) for which we did not 
obtain a valid redshift. If such a companion exists, we as- 
sign to it the same redshift as the primary galaxy and require 
that the companion lies within the selection lines in Fig. 1. 
This identifies all pairs (a total of 3) which are potentially 
missed because of redshift incompleteness. Pairs where both 
members are missed by the spectroscopic survey, will share in 
the general 2SLAQ incompleteness, without a bias for incom- 
pleteness at small angular separations. We can assume that 
all these 3 'extra' pairs are real and treat them as a source of 
systematic error on our determination of the pair fraction and 
the merger rate. Figure 3 shows postage stamp images of the 
two dynamical pairs we find and of the three possible pairs. 

P02 also use a weight Wg which accounts for the local mag- 
nitude incompleteness around each galaxy, a geometric effect 
due to slit placement and limiting filters in the CNOC2 survey, 
a color term and a term that depends on the evolution of the 
galaxy luminosity function over the redshifts covered by the 
survey. This weight is specifi c to the methods employed by 
the CNOC2 survey (lYee et alj ri996). In our case, we have a 
very homogeneous sample, there is no geometric effect (other 
than the one corrected by wg), all galaxies have similar colors 
and the color selection is very efficient in selecting galaxies of 
the appropriate magnitude and redshift. Eve n in the CNOC2 
survey most of these weights are equal to 1 (lYee et alJI 19961) 
for most galaxies. For these reasons, we do not adopt this 
weight in our study. 

3.5. Calculation of pair fraction 

Following POO the number of close companions per galaxy 
in a flux-limited sample can be expressed as: 
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Fig. 3. — Images (from the SDSS) of pairs found in our survey. The two 
left-hand ones are the real ones, while the three images on the right hand 
show the 'possible' pairs 
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and the total companion luminosity: 



(4) 



(5) 



The sums run over the i — 1, ...,Ni primary galaxies: Nc- 
and Lc- are the number and summed luminosity of galaxies 
from the secondary sample that are dynamically close (as de- 
fined above) to the i*^' primary galaxy and are expressed as: 
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where the sums run over those j — 1, ...,N2 galaxies in the 
secondary sample that fulfill the criteria for being dynamically 
close to the i*'* galaxy in the primary sample. 

The weights to be applied to the secondary sample correct 
for spatial incompleteness (wl,wi) and the flux-limit of the 
survey {S{zj)) for each j*'' companion. We do not apply the 
spectroscopic incompleteness weight wg because we use the 
alternate method described above to calculate the contribution 
from missed close pairs. The weights are defined in the pre- 
vious subsections. Similar weights also need to be applied to 
the primary sample, to correct for spatial incompleteness and 
the mean density: 
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1/2 (i.e., the reciprocals of the 



weights applied to the secondary sample). No spectroscopic 
completeness weight is needed for the primary sample. 

3.6. Merger timescales 

In order to convert a pair fraction into a merger rate the 
timescale for the merger event needs to be estimated. Merger 
timescales are difficult to determine and depend on the poorly 
known details of the merger process and how the poten- 
tial of the individual galaxies reacts to the merger episode. 
The more commonly used timescales in previous work and 
semi-analytic models are based on dynamical friction argu- 
ments. The dynam ical friction timescale can be calculated as 
jPatton et al.ll2000i) : 

2.64X IQSr^ 
= Mb^A ^^^^ 

where r is initial physical pair separation in kpc, Vc is the 
circular velocity in km s~^, M is the mass and A is the 
Coulomb logarithm. Assuming r ~ 20h^l kpc, Vc = 260 
km s^^ and In A = 2 for equal mass mergers, we get a typical 
dynamical friction timescale between 0.1 and 0.3 Gyrs over 
the range of masses we sample. 

The calibration by iKitzbichler & White! (|2009|) for the 
merger timescale of pairs in N-body simulations with M > 
5 X 10^ Mq and a velocity separation of 300 km s^^ yields 0.9 
Gyr, while the estimated merger timescale from the N-body 
simulations of iBovlan-Kolchin et al.l (l2008h is ~ 0.8 Gyr for 
the mass ratios we consider. A comparison between merger 
timescale by dynamical friction and in N-body simulations by 
[Bovlan-Kolchin et al. (2008) shows that the dynamical fric- 
tion formula underestimates the merger timescale by a factor 
between 1.7 and 3.3 for mergers of mass ratio 1:3 to 1:10 
respectively, and leads to a 40% overestimate of the mass ac- 
cretion rate (mainly via minor mergers). 

4. DISCUSSION 
4.1. Results 

We can now apply the above procedure to the 2SLAQ sam- 
ple. Out of 7889 galaxies we find a total of one dynamically 
close pair (we find a second pair, but it lies within the 3" ex- 
clusion circle). This yields a pair fraction Nc = 0.047% and 
a luminosity accretion rate Lc = 2.5 x 10*^ Lq. For galaxies 
within z < 0.55 this is Nc — 0.041%. Normally, errors on 
these quantities can be estimated by bootstrap resampling or 
jack-knifing, but this is not feasible with a sample of only two 
objects. We then proceed to estimate an upper limit to the pair 
fraction and the merger rate. 

In order to derive an upper limit to the merger rate, we take 
the observed detections and use the binomial distribution. For 
large samples (N > 100) the error distribution is given by 
dBurgasser et al.ll2003h : 



(11) 



where £{, is the pair fraction, n is the number of pairs, N 
is the number of objects in the sample, is the upper la 
probability limit to the pair fraction and is the lower la 
probability limit to the pair fraction. The region between 
and corresponds to the 68% confidence interval for ei, for a 



6 



De Propris et al. 



Gaussian distribution. In order to obtain aconservative upper 
limit we decide to use both pairs that we actually find (even 
though one of these is not within the valid range of separa- 
tions). 

For n = 2 and iV = 7889 this yields = 0.094±0.066 %. 
The systematic error due to the three possible extra pairs is 
0.14 ± 0.08%. We can therefore quote a 5cr upper limit to the 
pair fraction of 0.44% (random) plus 0.54% (systematic), for 
a total of 1.0%. 

As a check, we estimate the pair fraction in each of the 
separate 2SLAQ survey patches: while most regions have no 
pairs, in one area we find a pair fraction of 0.30%, which 
is consistent with the upper limit we derived above. Fi- 
nally, we can also consider pairs with wider separation in 
both projected distance and velocity: rp < 50 h^^ kpc and 
Av < 1000 km s^^. This considerably increases the num- 
ber of contaminants (u nphysical pairs), as shown by POO and 
iDe Propris et al.l (l2007h . However, this pair fraction may pro- 
vide a further interesting constraint. For the 2SLAQ sample 
we fmdN, = 0.13%. 

Our upper limit to the pair fraction can be translated into a 
5(7 upper limit to the dry merger rate using the expression: 

Rmg = Nc n{z) 0.5prnerg T^] (12) 

where n{z) is the space density of galaxies in our sample, 
0.5 is a factor introduced to avoid double counting galaxies (1 
pair contains 2 galaxies), Pmerg is the probability that pairs 
will merge (assumed to be 1 here) and Tmg is the merger 
timescale, for which we take the shortest possible timescale 
set by dynamical friction. We find that a robust 5(t upper 
Hmit to the merger rate is: < 0.8 x 10~^ Mpc~'^ h^^ Gyr~^ 
(including the systematic contribution due to the three extra 
photometric pairs) for galaxies with —23 < Mj. < —21.5 
at 0.45 < z < 0.65 (a period of 1.4 Gyrs in the history of 
the Universe). A more realistic merger timescale and merger 
fraction may decrease this estimate by more than one order of 
magnitude. 

4.2. Comparison with Previous Measurements 

The low dry merger rate we measure here is i n good agree- 
ment with a n umber of other estimates: locallv. lMasiedi et aH 
(I2006L I2008h obtain an upper Hmit of < 1.7% per Gyr for 
the dry merger rate of SDSS LRGs with Mi < -22.75 and 
z < 0.36, a value which is consistent with ours, albeit for 
more luminous (massive) galaxies. This is al so similar to the 
dry merger rate of 0.8% measured bv ,Wen et al.l (120081) for 
SDSS LRGs with Mr < -21.5 and z < 0.12. The inte- 
grated dry merger rate for galaxies in the Red Sequence Sur- 
vey is ^ 6% since z = 0.8 over —25 < Mr < —20, while 
iBundv et ai] (l2009l) find very low to zero likely dry mergers 
in GOODS data at z < 1.2. 

Nevertheless there are some discrepant estimates in the lit- 
erature. Apparently, the most worrying is the value of 2.4% 
pe r Gyr for 0.19 < z < 0.55 LRGs from 2SLAQ derived 
bv lWakeetan(l20()8h using the small scale correlation func- 
tion. However, IWake et alJ f2008 ) use LRGs over the en- 
tire range of absolute luminosities in the 2SLAQ survey, and 
their sample therefore includes both minor mergers (lumi- 
nosity rati o greater than 1:4) and less luminous obje cts. A s 
shown by Patton & Atfield ' (200S) and de Ra vel et al.l (l2008f) . 
the merger rate increases significantly for minor mergers and 
less luminous galaxies. Therefore the larger value derived 



bv IWake et al.l (l2008h is not necessarily in disagreement with 
ours. 

iLin et alj (l2008 l) apply the same method as we do (dynami- 
cally close pairs) to galaxies in the DEEP2 survey and derive 
an integrated merger rate of 24% since z < 1.2 for galaxies 
with —21 < Mb < —19. Assuming a flat evolution of the 
merger rate dLin et al.ll2008h . this is equivalent to ~ 6% per 
Gyr. Here, the different luminosity ranges sampled, the band- 
pass difference, and the fact that dry mergers were selected by 
morphology, rather than by colors and spectral features as we 
do, ma y play a role in explaining the difference. 

Bell et al.l (l2006ah obtain an integrated dry merger rate of 
^ 80% since z = 0.8 for galaxies in the Combo-17 survey 
with Mb < —20.5, while apply ing the small sca le corre- 
lation function to the same data, IBell et al.l (l2006b) derive a 
merger rate of 4% per Gyr at 0.4 < z < 0.8. As noted by 
IWake et alJ (120081) . the space density of objects in these sur- 
veys is 20 times greater than ours, and therefore iBell et all 
(2006a.b) are sampling considerably less luminous objects 
than we do. Given the dependence of the merger rate on lumi- 
nosity (Patton & Atfield 2008; de Ravel et al. 2008) this does 
not mean that our resu lts are in disagreement. In addition, 
iKhochfar & SiUd ilOOS h show that the above sample includes 
both wet and dry mergers, unlike ours where we confirm by 
spectroscopy that the vast majority of the sample has LRG- 
type spectra (Roseboom et al. 2006). 

Finally, Whi te et al.l ([2007) measure an integrated dry 
merger rate of 30% for LRGs in the NDFWS survey between 
0.5 < z < 0.9. This is equivalent to a merger rate of 3.4% per 
Gyr, but needs to be corrected for the higher mean redshift and 
different space density than the 2SLAQ sample. Wake et al] 
(l2008l) estimate that this downwa rd revision i s about a factor 
of 10, which places the result by IWhite et al.l (i2007.) in good 
agreement with ours. 

In addition, there are several estimates of the merg er rate 
for aU galaxies using asymmetriesjfor the pEEP2 ( Lot z et all 
120081) and COSMOS surveys dKampczvket a l. 2007), as well 
as galaxies in the GOODS fields dLopez-S aniuan et al. 2009), 
and pairs in the COSMOS survey (iKaitaltepe et al. 2007). 
The measured merger rate is approximately 2-4% per Gyr 
for the entire population of L > L* galaxies, of which dry 
mergers are only a subset, which is in agreement with our 
measurement. 

Our result is also consistent with the more indirect merger 
fraction derived from the evolution of the red galaxy lumi- 
nosity function. Note that, in this case, 'growth' may take 
place via the transformation of blue galaxies into quiescent 
objects, moving on to the red sequence, but without (necessar- 
ily) any major merging. IWake et aTl ( l2006l) used the 2SLAQ 
data to show that the LRG luminosity function evolves pas- 
siv ely at z < 0.6, a result confirnied by several other stud- 
ies (Bundv et al. 20061: ICaputi et al.l 2006: Cimatti et al. 200d 
Scarlat a et ak 2007). For ^ AL* LRGs Bjown et aL QMA 
2008) and lCool etaP (l2008h show that these galaxies evolve 
essentially passively at z < 1. 

Some massive red mergers are observed in distant clus- 
ters and groups (e.g., van Dokkum et al. 1999, 2001). Us- 
ing the pair fraction, van Dokkum et al. (199 9) estimated that 
^ 50% of massive red galaxies in the cluster MS 1054-0321 
at z = 0.82 ha ve undergone a merger to the present epoch. 
iMcIntosh et alJ (2008) calculate that, for groups and clus- 
ters in the SDSS at z < 0.12, the merger rate for mas- 
sive red galaxies is 2-9% per Gyr. These are much higher 
merger rates than we measure, albeit in a much denser en- 
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vironment. Assuming that the merger rate scales with den- 
sity, and since clusters are about 100 times denser than the 
fiel d, the values de r ived b y Ivan Dokkum et all ([T999L l200Th 
and lMcIntosh et al] ( l2008l) are probably not fully inconsistent 
with our ^ 0.3% upper limit for field LRGs. 

4.3. Implications for Galaxy Formation 

At face value, our results are in considerable ten- 
sion with models of galaxy formation in ACDM cos- 
mologies where most massive red g alaxies grow via dry 
mergers ([K hochfar & Burkerl! '2003'; 'De Lucia eXsD 120061 : 
iKhochfar & Silk 2008). Khochfar & Silk (2008) use a semi- 
analytic model to model the quenching of star formation and a 
standard merger tree, to predict the dry merger rate for galax- 
ies above a characteristic mass (the mass for which star forma- 
tion is shut off) at z < 1. For galaxies with M > 6.3 x 10^° 
M0 IKhochfar & Silk! (l2008h predict a dry merger rate of 
6 X 10~^ Mpc~^ Gyr~^, and a typical merger fraction 
of about 10-20% per Gyr, almost independent of redshift. 

Our luminosity range of —23 < Mr < —21.5 translates 
into a mass range of 16 to 5 xl0^° Mq. The upper limit to 
the merger rate we derive (< 0.8 x 10^^ Mpc^^ h^^ Gyr^^, 
including the 'worst case scenario' for systematic errors and 
assumuing all pairs merge on the shortest possible timescale) 
is at least a factor of 5 below these predictions. It is there- 
fore unlikely that dry mergers at z < 0.7 are important in the 
buildup of the red sequen ce (iBundy et al.ll2007l : iGenel et aP 
I2008bl : lBundvetani2009h 

On the other hand, the mass fraction in < L* galaxies 
on the red sequence appears to grow by about 50% since 
z = 1 (lBelletalJl2004l:lBorchetal.ll2(M 
ICool et al.ll2008l) . IScarlata et all (l2007r also note a deficit of 
fainter early-type galaxies in their COSMOS data at z = 0.7. 
The main mode of growth of the red sequence at z < 1 may be 
via cessation of star f ormation a nd morpholo gical evolu tion in 
lower mass galaxies jBell et ai.u2004 ; Faber etai.ll2007h . 

We therefore favor a scenario where most massive galax- 
ies are formed quasi-monolithically at high redsh ift and major 
mergers are relati vely unimpo rtant at z < 1 (e .g., [ Bower et al] 
120061 : iNaab et al.H2007i - cf., IConselicd l2006 for an observa- 
tional perspective). This is consistent with recent CDM simu- 
lations, that downplay the role of major mergers in galaxy for- 
mation over the past half of the Hubble time and favor minor 
mergers and internal proc esses as drivers of galaxy evolution 
in the last ^ 8 Gyr (e.g ., ICattaneo et alll200 8': 'Guo & Whitg 
120081 : iParrv et al.ll2008b . However, even in these models the 



more massive LRGs should grow primarily by major dry 
mergers, a result that appears to be in so me conflict with the 
slow growth observed for these galaxies jBrown et alj|200'7l 
l2008l : ICool etaPIIOOSh . 

It is important to distinguish between the merger rate of 
dark matter halos and that of their galaxy tracers. Strictly 
speaking, theory can only predict the former, while the lat- 
ter depends, at least in part, on comple x details of baryonic 
physics and dynamical friction (cf. Ber rier et alj |2006). Al- 
though our results appear to favor a particular realization of 
ACDM models, interpretation must await more detailed and 
realistic simulations of the stellar components of galaxies and 
their fate during mergers and interactions. 
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